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ABSTRACT 

We simulate the cooling of the neutron star in the X-ray transient KS 1731-260 after 
the source returned to quiescence in 2001 ft-om a long (> 12.5 yr) outburst state. 
We show that the cooling can be explained assuming that the crust underwent deep 
heating during the outburst stage. In our best theoretical scenario the neutron star 
has no enhanced neutrino emission in the core, and its crust is thin, superfluid, and 
has the normal thermal conductivity. The thermal afterburst crust-core relaxation in 
the star may be not over. 
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1 INTRODUCTION 

KS 1731-260 is a neutron star X-ray transient whose 
. observational histo ry has been described recently by 
' ICackett et all (|2006l ). The source was discovered in the ac- 
, tive state in August 1989 by the Kvant orbital observatory; 
' subsequent ana lysis showed that it had also been active in 
October 1988 (|Sunvaev et al.l [19901 ). It remained a bright 
, X-ray source showing type I X-ray bursts for about 12.5 
• years. It is believed that this activity was powered by accre- 
tion onto the neutron star (through an accretion disk) from 
. its low-mass companion in a compac t binary. Adopting a 
' distance to the source of D = 7 kpc, ICackett et alT (|2006l ) 
report a characteristic 2-10 keV luminosity of KS 1731-260 
in the active state of ~ 10'^^ erg s~^. The source remained 
active till the beginning of 2001 and then returned to qui- 
escence. The last detection in the active state was made 
on January 21, 2001 with RXTE, but on February 7, 2001, 
RXT E already failed to de tect KS 1731-260 in the active 
state (|Wiinands et al.ll2001^ . 

The first detection of KS 1731-260 in quiescence was 
made bv lWiinands et abl (|200ll ) with Chandra on March 27, 
2001. For D =7 kpc, the 0.5-10 keV luminosity was ~ 10^^ 
erg s~^, three orders of magnitude lower than in the active 
state. The radiation spectrum contains a component that 
can be interpreted as the thermal emission from the neu- 
tron star surface. Since then the source has been observed 
several t imes with Chandra a nd XMM-Newton as summa- 
rized by ICackett et al.l (|2006l ). Its X-ray light curve faded 
over time scale ~ 2 years showing a trend to flattening (with 
the residual luminosity of ~ 2 x 10^^ erg s~^). 

According to observations, the accretion in quiescent 
states of X-ray transients is stopped or strongly suppressed. 



The natur e of quiescent X-ray emission is a subject of de- 
bates (see. ICackett et al]|2006l and references therein, for a 
list of possible hypotheses). Here, we focus on the hypoth- 
esis of deep crustal heating of neu tron stars proposed by 
iBrown. Bildsten fc Rutledgel (|l998l ). It states that, when a 
neutron star accretes, its crust is heated by nuclear trans- 
formations (mainly by beta captures and pycnonuclear re- 
actions) in the accreted matter sinking within the crust un- 
der the weight of newly accreted material. The star remains 
sufficiently warm after an accretion episode, producing qui- 
escent surface emission. The sequence of nuclear transfor- 
mations and associated energy generation rates were cal- 
culated bv iHaensel fc Zd unik (199Qi) assuming that the ac- 
creted matter burns to ^®Fe in the neutron star surface layers 
so that, initially, before sinking within the deep cru st, the 
matter is composed of ^^Fe. Later lSchatz et al] (|200ll ) calcu- 
lated explosive nucleosynthesis in the neutron star surface 
layers and showed that the explosive b urning can proceed 
to mu ch heavier elements. Accordingly, iHaensel &: ZdunikI 
(|2003| ) proposed new deep crustal heating scenarios (start- 
ing with heavier elements, particularly, w ith ^"^Pd). In all 
the cases 'Hacnscl fc ZdunikI l|l99d . |2003| ) obtained similar 
deep crustal energy release, ~ 1 — 1.5 MeV per one accreted 
nucleon, sufficient to powe r quiescent therma l emission in 
X-ray transients. Recently iGupta et al.l (|2007h have recon- 
sidered the heating starting with multicomponent matter 
(ashes of explosive burning in the surface layers). They have 
shown that the heating of the deep outer crust can be higher 
because beta captures can produce daughter nuclei in ex- 
cited states; their deexcitation can generate extra heat. 

The onset of the quiescent state of KS 1731-260 was 
recognized as an outstanding phenomenon from the very 
beginning. The majority of other X-ray transients undergo 
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short accretion episodes (days to months) in which the deep 
crustal heating cannot break the crust-core thermal coupling 
and make the crust much hotter than the stellar core. How- 
ever, it is p ossible in KS 1731-26 because of the long accre- 
tion stage l|Rutledge et al.ll2002l ). Therefore, observations of 
its quiescent thermal emission can help to understand how 
the crustal heat spreads over the entire star, that is useful 
for exploring the neutron star structure. 

The first modelling of the KS 1731-260 cooling was 
made by iRutledge etahl l|2002l ) soon after the quiescence 
onset. The authors based o n previous simulations by 
lUshomirskv fc Rutledgj |20o3) of the crust-core relaxation 
in a neutron star with a heated crust. iRutledge et~aLl (|2002l ) 
proposed several c ooling scenarios bas e d on t he deep crustal 
heating model of iHaensel fc Zduni^ (Il990l ) and different 
crust and core microphysics. They predicted that the neu- 
tron star can reach the crust-core relaxation and associated 
flatte ning of the qu i escen t soft X-ray light curve in 1-30 
years. ICackett et al.l (|2006l ) have compared t he new observa- 
tions of KS 1731-260 with the predictions of lRutledge et al.l 
l|200^ and conclude that the star should have high thermal 
conductivity in the crust and enhanced neutrino emission in 
the core. 

Here we present new cooling calculations and discuss 
their consistency with the observations of KS 1731-260. 



2 COOLING MODEL AND PHYSICS INPUT 

Our cooling simulati o ns a re similar to those o f 
lUshomirskv fc Rutledgj l|200lh and IRutledge et all l|2002h . 
We assume that the neutron star crust underwent deep 
crustal heating during the long accretion stage. We employ 
the m odel of deep crustal heating of iHaensel fc ZdunikI 
l| 19901 ). but modify the energy release due to sequences of 
pairs of beta-captures in the crust. Specifically, we assume 
that daughter nuclei after a primary beta capture are 
produced in excited states and deexcite before a secondary 
beta capture, heating thus the matter (instead of wasting 
extra energy into neutrino emission). In this way the 
distribution of heating sources remains the same as in 
IHaensel fc ZdunikI (|l990D but the sources in the outer crust 
become stronger, resembling those obtained bv I Gupta et al.l 
||2007^. The sourc e posit ions and strengths, calculated by 
IHaensel fc ZdunikI l|2007h . are shown in Fig. [T] The overall 
energy release is 1.9 MeV per accreted nucleon. 

To simulate the neutro n star cooling we use our gen- 
eral r elativistic cooling code (|Gnedin. Yakovlev fc PotekhinI 
I2OOII ). It solves the thermal diffusion problem within the 
star (at densities p > ph) and uses a predetermined quasi- 
statio nary relation Ts — Tb (jPotekhin. Chabrier fc YakovlevI 
1 19971 ) between the effective surface temperature Ts and 
the temperature Tb at the base (p = pb) of a thin heat- 
blanketing envelope (p ^ pb). Now we shift pb from previ- 
ously used values ~ 10^" — 10^^ g cm~^ to pb = 10* g cm""'. 
This allows us to put all heat sources into the region of 
p > Pb and to reduce the time of heat propagation through 
the blanketing layer from ~ 1 yr to ~ 1 d (enabling the 
code to trace short-term - 1 day - surface temperature vari- 
ations) . 

To explore the sensitivity of calculations to the crust 
physics, we employ two models of the neutron star crust. 



composed of ground-state (GS) or accreted (A) matter. 
The g round-state crust (e.g.. IHaensel. Potekhin fc YakovlevI 
l2007h has been u sed in our pr evious simulations. The model 
of accreted crust (|Haensel fc Zdunik.199 01 is consistent with 
the adopted model of deep crustal heating. The accreted 
crust is composed of lighter nuclei with lower atomic num- 
bers. Deep in the inner crust, at p > lO^"' g cm~^, composi- 
tion is similar to the ground-state one, with >80% of nucle- 
ons constituting a neutron gas (Hacnsc l fc Zdunik. 1990 ). 

We employ the electron ther mal conductivity i n the 
crust, K, limited by electron-ion fenedin et al.l I2OOII ) and 
electron-electron (|Shternin fc Yako vlev 200^ scattering. It 
will be called normal. We will als o use the mod el electron 
thermal conductivity proposed by iBrownl 1 200^. It corre- 
sponds to an amorphous crust fe.g.. Ijones I2OO4I ) and will 
be called low. Actually, it gives the lowest limit on k in the 
crust. Several model thermal conductivities as functions of 
density in the crust for two values of temperature (T = 10* 
and 10^ K) are plotted in Fig. [T] 

In the inner crust, we take into account the effects of 
neutron superfluidity on the heat capacity of free neutrons 
(e.g., Yak ovlev, Levenfish fc Shibanov 1999i'). A representa- 
tive set of models for superfluid neutron gaps in the in- 
ner crust, which determine s uperfiuid critical temperatur e 
profiles Tc(p), is collected bv iLombardo fc Schulzd l|200ll ). 
The collection includes a well defined gap provided by the 
pure BCS theory of singlet-state neutron pairing (with a 
maximum of Tc ~ 2 x 10^" K within the crust) and a 
number of gaps calculated using various neutron polariza- 
tion models (with the maxima of Tc approximately three 
times lower). BCS superfluidity very strongly suppresses 
the neutron heat capacity in the inner crust; this super- 
fluidity will be called strong. The effects of other superfiuid 
models are weaker and more or less similar. For illustra- 
tio n of the latter effects, we will use t he model proposed 
bv lWambach. Ainsworth fc Pinei l|l993l ): such superfiuidity 
will be called moderate. We calculate the ne utrino emission 
in the crust and in the core according to lYakovlev et al.l 
(2OOII ). In our cooling models the neutron star stays not too 
hot, so that crustal neutrino emission (including that due to 
Cooper pairing of neutrons) is insignificant. 

In the neutron star core, we use an equation of 
state of dense matter (containing nucleons, electrons, and 
muons) constructed bv lAkmal. Pandharipande fc Ravenhalll 
()199& ) (their model Argonne V18-|-5«-|-UIX*). Specifi- 
cally, we adopt its convenient para metrization proposed by 
iHeiselberg fc Hiorth-JeiisenI (|l999l ) and described as APR I 
bv lGusakov et al.l (12005^ In this case, the maximum grav- 
itational mass of stable neutron stars is Mmax = 1.923M0 
and the direct Urea process of powerful neutrino emission 
opens at M > 1.828Mq. We will mainly use two neutron 
star models, with masses Ai = 1.6 and 1.4 Mq, where di- 
rect Urea process is forbidden; both stars demonstrate slow 
neutrino cooling via the modified Urea process. The 1.4 
star has the central density pc — 9.4 x 10^* g cm~^, the 
circumferential radius R = 12.14 km, and the crust thick- 
ness AR — R — Rcoia = 1.16 km (where Rcoie is the 
core radius corresponding to p = 1.5 x 10^* g cm~^). The 
1.6 Mq star is more compact, with thinner crust, and has 
Pc = 1.16 X 10^^ g cm"^, R = 11.88 km, and AR = 890 m. 

The thermal c onductivity of the neutron s tar co re is 
described following iBaiko. Haensel fc YakovlevI l|200lh and 
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Figure 1. Density dependence of the electron thermal conductiv- 
ity K (left vertical scale) in the neutron star crust with accreted 
(A) or ground-state (GS) matter at two temperatures (logT[K]=7 
and 8, numbers next to curves). The thin lower curve is for the 
model of low k while other curves are for normal k. Vertical bars 
show positions and power (right vertical scale) of the heat sources. 
Initial layer is assumed to consist of ^^Fe, as in lHaensel fc Zdunitd 
b ut neutrino losses in electron captures are suppressed, 
following [Gupta et al.l ||2007| '). 



IShternin fc YakovlevI l|200'i1 ). For simplicity, the effects of i 
cleon superfluidity in the core are neglected. 



3 RESULTS AND DISCUSSION 

We have calculated (Fig.[2J a number of cooling curves which 
give the effective surface temperatures , as detected by 
a distant observer, versus time t; t = Q refers to February 
1, 2001, the date near which KS 1731-260 turned in qui- 
escence. We compa re the curves with s even observational 
points presented bv lCackett et al.l (|2006l ): the values of 
were inferred from the observed X-ray spectra (employ- 
ing non-magnetic neutron star hydrogen atmosphere models 
from the Xspec database and assuming D = 7 kpc, i? = 10 
km, and M — IAMq). We doubled the reported la obser- 
vational error bars to enlarge statistical significance, that 
would make our analysis more realistic. 

To start any cooling calculation, we have taken a neu- 
tron star model with the thermally relaxed interior and 
some initial surface temperature . Then we switch on 
deep crustal heating produced by a constant mass accretion 
rate Ad over 12.5 years. In that period a certain amount of 
heat Etot is deposited into the crust. The crust is heated 
and its thermal balance with the thermally inertial core is 
violated. Then we switch off accretion (deep crustal heat- 
ing) and the crust cools down regaining thermal equilibrium 
with the core. Some (typically small) part of Etot diffuses 
to the surface and radiates away via thermal surface emis- 
sion. The rest is carried by thermal conduction to the core. 
The core temperature stays almost unchanged because of 
the high core thermal conductivity and heat capacity. The 
crust-core thermal relaxation takes 1-100 years, depending 
on the neutron star model. After this relaxation is over, the 
surface temperature nearly reaches its initial value . The 
star cools down further with isothermal interior over typical 



cooling time scales 1-10 kyr until the next accretion episode. 
The extra heat deposited to the core is mainly emitted over 
those long cooling time scales via core neutrino emission. 

Our cooling curves in Fig. [2] are calculated for differ- 
ent neutron star masses, microphysics in the crust, mass 
accretion rates (and -Etot), and (as shown in the fig- 
ure and Table [l| . Our aim is to explain the observed tem- 
poral evolution T^{t) of KS 1731-260 in the quiescent 
state. A successful explanation should also be consistent 
with the observational constraint on the mass accretion rate, 
M < 5 X IQ-^Mp, yr~^ (for D = 8 kpc, see Table 3 in 
lYakovlev. Levenfish fc Haense]|[2003l ), which translates into 
Etot ^ 2.4 X 10** erg for the adopted deep heating model. 
Table [1] shows that all presented cooling models roughly sat- 
isfy this requirement. 

Fig. [2^ refers to the 1.6 Mq neutron star model, while 
Fig. [2]3 is for the IAMq star. All curves in Figs. [2^ and 
b are calculated assuming the initial surface temperature 
to be = 0.8 MK (so that the internal temperature is 
~ 8 X 10^ K, as in a cooling isolated neutron star which 
is ~ 10^ years old). This is a typical surface temperature 
of the neutron star provided by the last three observational 
point s. Thus, in Figs. [2^ and b we (following ICackett et al.l 
I2OO6I ) tacitly assume that the crust-core equilibrium is re- 
established in two years after the quiescence onset. In all 
curves in Figs. [2^ and b, but in curve 6, Etot has been cho- 
sen in such a way for the surface temperature at the first 
quiescent observation to be consistent with data. 

Curve 1 in Fig. [2^ seems to be the best. It corresponds 
to the accreted crust with the normal thermal conductivity 
and moderate neutron superfiuidity. It naturally explains 
the thermal relaxation of KS 1731-260 with the standard 
physics input. The maximum internal temperature raise to 
r ~ 4 X 10* K takes place at t — near the bound- 
ary between the outer and the inner crust. The core-crust 
relaxation takes ~ 2 years. The star would need ~ 10^ 
years to emit all the heat pumped into the core during the 
outburst and reach the same thermal state as before the 
outburst. This is in good agreement with the estimate of 
iRutledge et"al] (|2002l ) for a similar cooling model. 

Using the same physics as for curve 1 but the ground- 
state crust (with lower conductivity) we obtain slower re- 
laxation (curve 3). It is acceptable but less consistent with 
the observations; it requires lower Etot because it is easier 
to heat the crust with smaller thermal conductivity. Taking 
the latter cooling model 3 and neglecting superfluidity in 
the inner crust we obtain curve 2. A non-superfluid crust 
has larger (neutron) heat capacity which noticeably delays 
the thermal relaxation making it much less consistent with 
the data. Returning to our best model 1 but assuming strong 
superfluidity in the crust, we stronger suppress the heat ca- 
pacity of neutrons and obtain curve 4; it shows faster and 
quite acceptable relaxation. The effects of strong and moder- 
ate superfluidity are actually very close, although the pres- 
ence of superfluidity greatly improves the agreement with 
the data. Now if we return to model 1 but assume low ther- 
mal conductivity, we come to much longer crust-core relax- 
ation (over several hundred years, curve 5). It is inconsistent 
with the obse r vation s, in agreement with the conclusion of 
ICackett et all l|2006l ). Finally, if we take the best model 1 but 
assume the same (lower) mass accretion rate as in model 2, 
we get curve 6. Therefore, the latter mass accretion rate. 
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Table 1. Cooling curves in Fig. 2 



Curve Tjfg' Crust Conduction Superfluid -Etot 

MK model in crust in crust 10*'* erg 



la 0.8 A normal moderate 2.6 

2a 0.8 GS normal none 1.9 

3a 0.8 GS normal moderate 1.8 

4a 0.8 A normal strong 2.6 

5a 0.8 A low moderate 0.6 

6a 0.8 A normal moderate 1.9 



lb 0.8 A normal moderate 2.3 

2b 0.8 GS normal none 1.7 

3b 0.8 GS normal moderate 1.5 



Ic 0.67 GS normal none 2.4 

2c 0.63 GS normal none 2.4 



being used for the microphysics of model 1, is insufficient to 
explain high values of in the beginning of the quiescent 
state. 

Curves 1-3 in Fig.[2j3 are analogous to curves 1-3 in Fig. 
[2^, but are calculated for a less massive star, with thicker 
crust. The thicker crust produces longer thermal relaxation, 
less consistent with the observations. 

We have also performed many other cooling calcula- 
tions varying physics input. In particular, we have varied 
the distribution of heat sources within the crust. We have 
obtained that it is much easier to explain the observations 
by placing the sources into the outer crust. These models 
naturally give short thermal relaxation and efficient heat- 
ing of the surface. In contrast, were all sources located in 
the deep inner crust, the star would show longer thermal 
relaxation and one would need too much energy to heat the 
surface because the heat would be pumped into the core. 
In connection to this, the improved model of deep crustal 
heating used here, where the heat release in the crust is en- 
hance d compared to the original model of lHaensel fc ZdunikI 
(|l990l ) (due to sw itching-off neutrino losses associated with 
electron captures. iGupta et al-llioOTl ). is more favorable for 
explaining the observations. 

In addition, we have artificially varied the thermal con- 
ductivity in different places of the crust and found high sen- 
sitivity of the cooling curves to these variations. The conduc- 
tivity strongly affects both, the thermal relaxation time and 
the efficiency of surface heating. Taking the conductivity a 
few times lower than the normal conductivity of accreted or 
ground-state crust produces too long crust-core relaxation 
which disagrees with the data. 

Furthermore, we have taken different neutron star mod- 
els (different equations of state in the core, and different 
masses). In particular, we have used the models of massive 
neutron stars whose core neutrino emission is strongly en- 
hanced by the nucleonic direct Urea process (e.g., 1.9 Mq 
model for the equation of state employed in Fig. [5]). We have 
found that we need unrealistically intense crustal heating 
(too high Etot) to explain the high observed values of T^{t) 
in the beginning of quiescence. Moreover, such a star has too 
short global cooling time scale (years to decades), compara- 
ble to the crust-core relaxation time. The crust-core relax- 



ation becomes coupled to the global thermal relaxation; the 
cooling curves do not show the observed flattening at f > 2 
years. Hence, we cannot reconcile theory with observations 
if the neutrino emission of the star is enhanced by the direct 
Urea process. Nevertheless, we think that it may be pos- 
sible to explain the observations if the neutrino emission is 
enhanced by a less efficient mechanism (e.g., by pion or kaon 
condensation in the stellar core) or if the direct Urea process 
operates but is strongly suppressed by nucleon superfluid- 
ity (e .g., l^kovlev fc Pethickll2004 IPage. Geppert fc Webeil 
l2006l ). 

Finally, we remark that the thermal crust-core relax- 
ation in KS 1731-260 may be still not over. This is illustrated 
in Fig. [2j:, where we present two new cooling curves for our 
1.6 Mq and 1.4 M© neutron star models. They are calcu- 
lated without imposing the constraint that = 0.8 MK. 
We have intentionally taken the physics input (ground-state, 
non-superfluid crust with normal conductivity) which gives 
too long thermal relaxation to explain the data for the sce- 
narios in Figs. [2^ and b. Now we take lower and reach 
consistency with the current observations (and get a rather 
low crustal heat release Etot). We see that the crust-core re- 
laxation in KS 1731-260 can really last longer than 2 years, 
and this possibility widens the class of cooling models con- 
sistent with the data. It will hopefully be checked in future 
observations of KS 1731-260. 



4 CONCLUSIONS 

We have simulated the cooling of the neutron star in the 
quiescent state of KS 1731-260 employing the model of deep 
crustal heating of the star in the outbu rst state. We have 
used the model of deep crustal heating l|Haensel fc Zduiii^ 
1990) up dated by switching-off neutrino losses in the crust 
(|Guptaeral]|2007l '). Our main conclusions are: 

(i) One can explain current observations of KS 1731-260 
using a model of deep crustal heating and a standard micro- 
physics of the neutron star. 

(ii) If the crust-core thermal relaxation in the neutron 
star is reached in ~ 2 years, the most successful cooling 
model implies the model of accreted crust with normal ther- 
mal conductivity and neutron superfluidity; the neutron star 
should be sufficiently massive (to have a thinner crust), but 
the neutrino emission in its core cannot be too high (e.g., 
it can be provided by the modified Urea process). All these 
factors shorten the crust-core thermal relaxation. 

(iii) The model of low thermal conductivity (amorphous 
crust) gives too long crust-core relaxation, inconsistent with 
the data. 

(iv) The enhanced neutrino cooling via the direct Urea 
process in the neutron star core gives too fast cooling of the 
entire star and requires too intense crustal heating, incon- 
sistent with the data. 

(v) The crust-core thermal relaxation can be not reached 
yet. If so, the data can be explained by a wider class of 
neutron star models. 

We stress that the thermal crust-core relaxation of the 
neutron star in KS 1731-260 is much more sensitive to the 
physics of the crust than the core. We employed the mod- 
els of non-superfiuid core just for simplicity. Core superflu- 
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Figure 2. Theoretical cooling curves for (a) M = 1.6 Mq and (b) 1.4 Mq neutron stars, and (c) for stars with both M compared with 
observations. The curves are explained in Table [T] and in text. 



idity can change the core heat capacity and neutrino lu- 
minosity, but the principal conclusions will be the same. 
Our calculations are not entirely self-consistent. For in- 
stance j_tlM_suTfa££tein£CT was inferred from observa- 
tions ( Cackett et al.l [2OO6I ) . assuming neutron star masses 
and radii different from those used in our cooling models. 
This inconsistency cannot affect our main conclusions, but 
it would be desirable to infer for our neutron star mod- 
els. The ther mal relaxation in the quiescent state has been 
observed also (ICackett et al.ll2006l ) for another neutron star 
X-ray transient, MXB 1659-29. We hope to analyse these 
data in the next publication. 



5 ACKNOWLEDGMENTS 

We are grateful to A. I. Chugunov, O. Y. Gnedin, K. P. 
Levenfish and Yu. A. Shibanov for useful discussions, and 
to the referee, Ulrich Geppert, for valuable remarks. The 
work was partly supported by the Russian Foundation for 
Basic Research (grants 05-02-16245, 05-02-22003), by the 
Federal Agency for Science and Innovations (grant NSh 
9879.2006.2), and by the Polish Ministry of Science and 
Higher Education (grant N20300632/0450). One of the au- 
thors (P.S.) acknowledges support of the Dynasty Founda- 
tion and perfect conditions of the Nicolaus Copernicus As- 
tronomical Center in Warsaw, where this work was partly 
performed. P.H. and A. P. thank the Institute for Nuclear 
Theory at the University of Washington for hospitality and 
the U.S. Department of Energy for partial support during 
the completion of this work. 



References 

Akmal A., Pandharipande V. R., Ravenhall D. G.,1998, 

Phys. Rev. C, 58, 1804 
Baiko D. A., Haensel P., Yakovlev D. G., 2001, A&A, 374, 

151 

Brown E. F., 2000, ApJ, 531, 988 

Brown E. F., Bildsten L., Rutledge R. E., 1998, ApJ, 504, 
L95 



Cackett E. M., Wijnands R., Linares M., Miller J. M., 
Homan J., Lewin W. H. G., 2006, MNRAS, 372, 479 

Gnedin O. Y., Yakovlev D. G., Potekhin A. Y., 2001, MN- 
RAS, 324, 725 

Gupta S., Brown E. F., Schatz H., MoUer P., Kratz K.-L., 

2007, ApJ, 662, 1188 
Gusakov M. E., Kaminker A. D., Yakovlev D. G., Gnedin 

O. Y., 2005, MNRAS, 363, 555 
Haensel P., Zdunik J. L., 1990, A&A, 227, 431 
Haensel P., Zdunik J. L., 2003, A&A, 404, L33 
Haensel P., Zdunik J. L., 2007, A&A, in preparation 
Haensel P., Potekhin A. Y., Yakovlev D. G., 2007, Neutron 
Stars. 1. Equation of State and Structure, Springer, New- 
York 

Heiselberg H., Hjorth- Jensen M., 1999, ApJ, 525, L45 
Jones P. B., 2004, MNRAS, 351, 956 

Lombardo U., Schulze H.-J., 2001, Lecture Notes Phys, 578, 
30 

Page D., Geppert U., Weber F., 2006, Nucl. Phys. A, 777, 
497 

Potekhin A. Y., Chabrier G., Yakovlev D. G., 1997, A&A, 
323, 415 

Rutledge R. E., Bildsten L., Brown E. F., Pavlov G. G., 

Zavlin V. E., Ushomirsky G., 2002, ApJ, 580, 413 
Schatz H. et al., 2001, Phys. Rev. Lett., 86, 3471 
Shternin P. S., Yakovlev D. G., 2006, Phys. Rev. D, 74, 
043004. 

Shternin P. S., Yakovlev D. G., 2007, Phys. Rev. D, 75, 
103004 

Sunyaev R. A. et al., 1990, SvA Lett., 16, 59 
Ushomirsky G., Rutledge R. E., 2001, MNRAS, 325, 1157 
Wambach J., Ainsworth T. L., Pines D., 1993, Nucl. Phys. 
A, 555, 128 

Wijnands R., Miller J. M., Markwardt C, Lewin W.H.G., 

van der Klis M., 2001, ApJ, 560, L159 
Yakovlev D. G., Pethick C. J., 2004, ARA&A, 42, 169 
Yakovlev D. G., Levenfish K. P., Haensel P., 2003, A&A, 

407, 265 

Yakovlev D. G., Levenfish K. P., Shibanov Yu. A., 1999, 

Physics-Uspekhi, 42, 737 
Yakovlev D. G., Kaminker A. D., Haensel P., Gnedin O. Y., 

2001, Physics Reports, 354, 1 



